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1. Introduction 

The dependence of critical exponents of statistical mechanical systems on the physical space 
dimension is one of the most fundamental features of critical phenomena [1] . Continuous changing 
of the space dimension d has been the basic principle for constructing the famous Wilson-Fisher 
epsilon expansion [2] with e = 4 — d. However, the Renormalization Group (RG) [3], which provides 
the fundamental theoretical basis for such calculations, is nonperturbative in its nature. It does 
not require that the critical exponents or any other universal quantities must be expanded in some 
space deviation, coupling constant, inverse number of order-parameter components, etc. 

While investigating the critical phenomena, the RG allows, in principle, to work in any setting 
of interest provided that we are smart enough to find out the appropriate analytical or numerical 
tools. The e-expansion has shown itself as a very useful tool in studying the qualitative features of 
systems in critical state even in its lowest-order approximations (see, e.g., [4]). Pushed to higher 
orders, it allowed to produce very accurate numerical extrapolations for critical exponents of three- 
dimensional TV-vector models [5] . An alternative calculational scheme used with comparable success 
is the so-called (^-expansion within the massive field theory in fixed dimension put forward by Parisi 
[6]. This approach has been mainly applied directly in three dimensions to systems of different 
complexity [7-12] (for a recent review see [13]). In [7] it was used for calculating the exponents of 
cj) 4 models in two dimensions, and in [14] — applied to pure and disordered Ising systems in general 
dimensions 2 < d < 4. 

A natural access to non-integer dimensions provides also the well-known large-AT expansion 
[15-18]. In its general scope, it gives explicit expressions for critical exponents as functions of d. 
These are valid in the whole range between the lower and upper critical dimensionalities, and can 
be handled analytically. The large-A^ expansion is capable to yield information on dimensional 
dependencies that are hardly accessible by other means, for example, from the epsilon expansion. 
Unfortunately, it is very hard to obtain such results in higher orders in l/N, while short series 
expansions in l/N usually fail to give very accurate numerical estimates for relatively small values 
of N, say N = 3. The convergence of truncated l/N expansions has been analyzed on the basis of 
the field-theoretical approach at d — 3 in [19]. 

The knowledge of dimensional dependencies of critical exponents, or eventually also universal 
relations of critical amplitudes are of great theoretical and practical interest. The information 
of this kind broadens out our fundamental knowledge about these main characteristics of the 
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critical behavior. Moreover, the possibility to consider the results, which are not directly related to 
dimensional expansions, in the vicinity of the upper or lower critical dimensions provides us with 
useful checks of their correctness or with some new relevant insights. 

Let us give several examples. At the beginning of 1980ies Newman and Riedel [20] used the 
exact Wegner-Houghton [21] RG equation and the scaling field method in a study of the pure and 
dilute Ising models (IM) in general dimensions from the interval 2.8 < d < 4. 

Pinn et. al. [22] studied the RG in the hierarchical model in 2 < d < 4. Some of their results 
have been given as tables illustrating the dimensional dependencies of critical indices calculated 
for some sets of discrete non-integer values of d. A contact with the e-expansion has also been 
established. 

Recently, Ballhausen et. al. [23] calculated dimensional dependencies of IM critical exponents 
for 1 < d < 4 by using the first-order derivative expansion of the exact RG equation for the effective 
average action (for recent reviews see [24,25]). Their results for v(d) and 77(d) are very similar to 
that obtained in [26] and [14]. 

In a very recent paper [27], O'Dwyer and Osborn analyze the Polchinski version [28] of the 
exact RG equations [21,3,24] and its possible truncations and derivative expansions. This is done 
both for non-integer dimensions between two and four, and in the epsilon expansion. Such com- 
bined approach allowed to better understand the local potential approximation [29] , the derivative 
expansions used in treatments of the Polchinski equation, and moreover, to suggest an alternative 
improved truncation of this equation. 

Another story happened in the study of the critical behavior of O(N) x 0{M) spin models 
[12,30-33]. Apparently, certain controversies in the discussion of these references could be avoided 
if there would be more well-established information related to the dimensional dependencies of RG 
functions, fixed points, and critical exponents. 

Finally, let us note that there is a big continuous interest in the high-energy physics (HEP) 
literature in calculating explicit closed-form results for Feynman integrals in general dimensions 
[34-45]. 

We believe that the results displaying the explicit dependencies of RG functions and physical 
observables on the space dimension d are of fundamental importance. The aim of the present paper 
is to do a little step in this direction. We return to the article [14] published some time ago with 
Holovatch and calculate the closed-form expressions for the two-loop RG functions of the massive 
field theory, treated in that work numerically, in general dimensions 1 < d < 4. 



2. Two-loop RG functions and Feynman integrals 



Let us consider the usual massive [46-49,6] 0(A^)-symmetric (f> 4 theory in general dimensions 
d < 4. To second order in renormalized coupling constant u, the corresponding (3 and 7 RG 
functions can be written as [49,50,14] 

(3(u) = -(4 - d)u{l -u+ [(57V + 22)/(d) + (N + 2)j(d)] } + 0( U 4 ), 



7 *(u) = -4(4-d) 
%2(u) = (4-d) 



N + 2 



{N + 8) 2 
iV + 2 



(N + 8) 2 
j(d) u 2 + 0(u 3 ), 



N + 8 



u- 12 



N + 2 

(TV + 8)2 



f(d)u 2 



+ 0(u 3 ). 



At the fixed point u = u* , which is determined by the zero of the j3 function f3(u) , the function 7^ 
gives the value of the Fisher exponent 77, while 7^2 (u*) leads to the combination 2 — v^ 1 — 77 where 
v is the correlation length critical exponent. 

The d-dependent functions f(d) = i(d) — 1/2 and j(d) are defined by the combinations of 
Feynman diagrams taken at zero external momenta, 



i(d) 



= X^/XZX 2 and M = A -Q- |^ /X ZX 



(1) 
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The massive two- loop integrals appearing in i(d) and j(d) are given by 

, 2 „ f f d d h d d k 2 1 

{m ' ' J J (2^) 2d (k 2 +m 2 ) 2 (k 2 +m 2 )[(k 1 +k 2 ) 2 +m 2 } ' W 

2 2 d) _ f f d d ki d d k 2 1 

J(m ,q;d)-J J (27r)M (fc2 + m2) (fc 2 + m2) [(fei + fc2 + g) 2 + m 2] ■ ( 3 ) 

The one-loop integral which appears in (1) and is traditionally used [51,52,7,50] for normalization 
of the coupling constant is given by 

D ^- d) " / wh? = "' J - 4 ^m - m ( 4) 

where T(z) is a usual Euler's gamma function (see e.g. [53]). This integral is absorbed, along with 
the symmetry factor (N + 8)/6 into the normalization of the dimensionless renormalized coupling 
constant u. With this normalization, the RG (3 function starts with — u + u 2 + 0(u 3 ) in the three- 
dimensional theory [51,7]. On the other hand, since the integral D{d) has a 1/e pole as d approaches 
4, in the epsilon expansion the fixed point value u* will be given by 1 + 0(e) as a consequence of 
such normalization. 

It is easy to see that the folowing relations hold true for the integrals I and J, 

I(m 2 :d) = -l-JL j(rn\ 0; d) = (m 2 ) d - 4 J(l, 0; d), (5) 

where the second equality immediately follows from the scaling property J(m 2 , q 2 ; d) = (m, 2 ) d ~ 3 J(l,q 2 /m 2 ; d) 
of the function J. 

In [14], using the Feynman parametrization (see, e.g., [54]), i(d) and j(d) have been reduced to 
the double integrals over Feynman parameters and written down as 

= r ( £ ) f 1 xdx C y e/2d y (f * 

() T 2 (e/2)J Q [ x (l-x)}^/ 2 J Q [x(l-x)(l-y) + y}^ {) 

T(e) r 1 dx r 1 v"Hl-y)dy 

J[ ' r 2 (e/2)J [x(l-x)}-^ 2 J [x(l-x)(l-y) + yY- [<) 

Here and further on, we parametrize the dimensional dependencies by the usual deviation from the 
upper critical dimension e = 4 — d. Generically, we do not assume that e is infinitesimally small. 

At the time of writing the paper [14], the following information was available. At d = 3, the 
integrals i and j could be easily evaluated analytically as i(3) = 2/3 and j(3) = —2/27. At d = 2, 
their numerical values are [52,50] i(2) = 0.78130241 ... and j(2) = -0.11463575 . . .. As d -» 4, 
both integrals 1(d) and J(d) have poles in s, namely 1(4 — s) <~ 1/e 2 and J (4 — e) ~ 1/e. In the 
functions i (4 — e) and j(4 — e) these poles are compensated by division through D 2 (4—e), so that 
their e expansions start with 1 [49] 



1 e 

i(d) = 

where 



M = j + j + 0( s ') and j(,i) = -|-£(2 + ;) + 0( e s ) (8) 



dx\- — } — r + n^'/r J> U ■ w 

(1 — x(l - x) [1 - x(l — x)\ z J 

This integral has not been calculated in [49] since it is "renormalization-dependent" and disappears 
from the results for the critical exponents [49]. Numerically, I = —1.5626048 . . .. For non-integer 
d, the functions f(d) and j(d) have been tabulated in [14] by using numerical computations in (6) 
and (7). The corresponding graphs have been plotted point by point for 2 < d < 4. The aim of 
the following section is to obtain analytically the explicit expressions for these functions in general 
dimensions d. 



x The similar formula given at p. 873 of [14] is not correct. 
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3. Explicit calculations in general dimensions 

Let us return to the integrals i(d) and j(d) given by (6) and (7) and consider the inner y integral 
from (6). Denoting for a while z = x(l — x) we write it as 



I f (z) = z- E f dyy e l 2 [l-y{z-l)/z\ 
Jo 



(10) 



We recognize that here we have to deal with an integral representation of the incomplete beta 
function B t (a,b) (see, e.g., [53], entry 6.6.1). But it is more convenient to express this function in 
terms of the Gauss hypergeometric function using the relation 6.6.8 of [53], 

B t (a,b) = a - 1 x a 2 F 1 (a,l-b;a+l;t). (11) 

Thus we obtain 

if(z) = z ~ £ 2*1 1 1 + £ A £ ; 2 + e / 2 ; (* - • (12) 

However, the argument of the resulting hypergeometric function is not good for further integration 
over x in the limits from to 1 since the factor 1/z is singular at x = and x = 1. The situation 
is substantially improved by using the three-term linear transformation formula 15.38 from [53] 
which reads 

2 F 1 ( a ,6;c;x) = ^| (1 - x)~ a 2 F 1 [a, c - b; a - b + 1; (1 - x)- 1 } 

+ a similar expression with a and b interchanged. (13) 
After some algebraic transformations we obtain for the inner y integral from (6) 

2 r -J^^^ 1 - £/2 (l-^r 1 - £/2 +2F 1 ( £ ,l;./2;,)" . 



W) = 



2-e 



2 m 



(14) 



An essential simplification occurred in the first term due to an appearance of a hypergeometric 
function with equal nominator and denominator parameters (see, e.g., [53], entry 15.1.8): 

2 F 1 (a,6;M = ]T^z fc = (l-z)- a (15) 

k>0 

where (a)& = T(a + k)/T(a) = a(a + 1) . . . (a + k — 1) is the Pochhammer symbol. 

Let us substitute the first term of (14) into the outer x integral remaining in (6). Recalling 
the short-hand notation z = x(l — x) we see that the factor z 1_£ / 2 is exactly cancelled by the 
analogous term in the denominator and, apart from the numerical constants, we remain with 

xdx _ f 1/2 dt _ 2 1+£ / £ 1 3 _1\ 

(I - X + X 2 )^/ 2 ~ Jo (3+i2)l+ e /2 ~ 2 l \ + 2' 2' 2' 3/" ( ' 

The last equality can be obtained by simply expanding the denominator of the t integral via (15) 
and integrating the resulting series expansion term by term. 

The second part of the square brackets in (14) can be handled in a similar fashion. A straight- 
forward calculation employing the Gauss series representation of 2-F1 yields 

f 1 xdx \ e , .1 r 2 (e/2) / 1 e 1\ 

Now we notice that the hypergeometric function with the argument —1/3 can be expressed in terms 
of a similar function with the argument 1/4 using the the linear transformation formula (see, e.g., 
[53], entry 15.3.4) 

2 F 1 (a,b;c;z) = (1 - z)- 1 2 F l (a, c - 6; c; -^-) . (18) 



4 



Two-loop RG functions 



Thus we get from (6) and (16) along with (18) and (17) our first explicit result 

1 



i(d) = 



£ . 1 

2' 4 



2^1 (l 



£31 

2' ' 2'4 



(19) 



As £ — > 0, the first hypergeometric function in square brackets reduces to 1, the second term 
vanishes due to the explicit factor s, and i{d) = 1/2 + O(s), as it should (see (8)). Despite of the 
overall factor 1/(2 — e), there is no singularity at d = 2 since the combination in square brackets 
vanishes at e = 2. In Appendix we derive the explicit result for i(2) in terms of the generalized 
hypergeometric function 3 F 2 . 

Successive application of Gauss' relations for contiguous functions 2F1 given by entries 7.3.1.16 
and 7.3.1.18 of [55] to the second hypergeometric function from (19) leads to an elegant represen- 
tation 



i(d) 



1 



iFi ( 



2Fl (l' 1; M) 



(20) 



Here an interesting symmetry of the nominator and denominator parameters in both functions 2-F1 
is observed, while they differ one from another by e/2. 

Let us consider the second function j(d). The calculation of the double integral (7) is somewhat 
more involved but the experience of the above calculation allows us to be brief. 

While considering the inner integration If(z) in (10) we have seen that it is useful to get its 
result in terms of a Gauss hypergeometric function of the argument z = x(l — x), and not (z—l)/z 
which is suggested directly by the integrals over y in (6)-(7) and (10). This can be achieved 
directly on the level of the integral representation. Indeed, changing the initial integration variable 
via y = l/(t + 1) we obtain for the inner integral Ij{z) from (7) 



- f 

Jo 



tdt 



( t + l)3-e/2 



(l + zt)- e = 



(4-£)(2-£) 



2 F 1 (s,2;e/2;z) 



(21) 



- z)- 2 -^ 2 [4 - e - 2(1 - e)z] . 



(4-£)(2-£) r(£) 

The two terms of Ij(z) are similar to that generated by the inner y integration in i(d) (see (14)). 
They can be handled in the same manner inside of the external integral over x in (7). According 
to this decomposition of Ij(z), we obtain two different contributions to the function j(d). We can 
write 



(4-£)(2-£) 



[h(d) + J 2 (d)} 



(22) 



where 



32(d) 



1 



3(l+e) 

2 l+e 
3^ 



2-Fl 



3+£ 1 



4(4- 



Fi 1 + £,1; 



5+£ 1> 



3(1 -e) 3(3 

/- N ( E 1 3 1\ „ „ / £ 1 3 1\ 

(1 - ^1 3; 2= -3) + C 1 + 2' 2 ; 2 5 "a) 

-|(2 + e) 2^1(2 



£ 1 3 
2' 2' 2' 



Similarly as before, using the linear transformation (18), we can reduce all above hypergeometric 
functions to that with the argument 1/4. Further, the repeated use of the contiguous functions' 
relations leaves us only with the same two 2F1 functions which appeared in i[d). Moreover, they 
build up just the combination defining that function. Finally we arrive at a simple expression of 
j(d) in terms of i(d): 

1 



(23) 



We recall that two equivalent explicit expressions for the function i(d) are given in (19) and (20). 
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4. Some results in two dimensions 



While the integrals i(d) and j(d) are smooth functions of d in the range of our interest, their 
calculation at d = 2 requires some care. In order to calculate i(2) from (19) or (20), we have to 
introduce a small deviation a from d = 2 and consider the limits a — > of these expressions. 
Details of the calculation can be found in appendix A. 

We obtain the explicit expression for i(d) at d = 2 in terms of the generalized hypergeometric 
function 3F2: 

/x 4/111331\ 
^ = ^ 3 M2'2<2 ; 2<2 ; l)- ^ 

Using the relation (23) it is easy to obtain the value of j(2): 

4 /l 1 1 3 3 1\ 2 
* 2) = - Vf 3F2 (2'2'2 ; 2'2 ; l) + 3' (25) 

Searching for some other representations of i(2) we have found the relation (details are given 
in appendix A) 

Here Cl2(7r/3) = 1.0149417. . . is the maximum value of the Clausen's function (for more informa- 
tion on this and related functions see [56]) 



Ch(9) = - f 
Jo 



dO In 



^- 2 \ = T, — - ( 27 ) 



Apparently, the identity (26) does not appear explicitly in the mathematical literature. It supple- 
ments the very similar formula 

3F2 V2'2'2'2'2' _ lJ " 10 ( 8) 

from mathematical tables [55], entry 7.4.6.1. 



5. What can we learn from the HEP literature? 

As well as in the condensed-matter physics, the calculation of multi-loop Feynman integrals is 
an important issue in the high-energy physics literature. Although the Feynman diagrams in four 
space-time dimensions are of main concern in the HEP theory, very often they are calculated by 
using the dimensional regularization in arbitrary space-time dimension d. It is recognized (see e.g. 
[34]) that working in general d dimensions is often much easier than directly at d = 4. 

Actually, the condensed-matter and high-energy physicists are often doing much the same work 
by pursuing physically different goals. This is, perhaps, most pronounced just in calculations of 
Feynman integrals. Unfortunately, it happens not very often that authors of these two communities 
pay enough attention to achievements gained in another one. The aim of the present section is to 
quote some HEP results related to calculations of the preceding sections. In doing this we show 
several useful implications of establishing such contacts. By reviewing briefly some HEP references 
we would like to draw attention of the cond-mat readers to them, as we feel that these are not well 
known and appreciated in the cond-mat literature. 

Some of the results obtained in [34] are directly related to ours. In particular, the authors of this 

reference calculate the two-point function J 3 (m 2 ,q 2 ) ~ — (^^) — . By using the Mellin-Barnes 

contour integral representation for massive propagators [57] they derive for J3(m 2 ,g 2 ) an infinite 
series over linear combinations of 3F2 functions 2 . 

2 Only quite recently an explicit expression for Jz{m?,q 2 ) has been found [40] in terms of Gauss and Appell 
hypergeometric functions. 
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The first two coefficients of the small-g 2 expansion of J^(l,q 2 ) are simply proportional to our 
functions i(d) and j(d). This can be seen by noticing (see (3)-(4)) that J(m 2 ,q 2 )/D 2 (m 2 ) = 
(e 2 /4) Jz(m 2 , q 2 ) and taking into account the relation (5). Thus we get, first, 



i(d) 



1 



2-e 3 



2-e 3 



3^ 



^(e) 3 /e 3.1x 

f2^72) + 2 (1_£)2Fl V2' 1 '2'4J 



(29) 



e/2. Analogous 



where we quote the equation (33) from [34] with the usual replacement e 
expressions can be found also in [37], eq. (4.38), and [58], eq. (A. 11). 

The last expression for i(d) is equivalent to that found in (19) and (20). A direct comparison 
between (29) and (19)-(20) can be done by applying the contiguous relations for hypergeometric 
functions with e/2 in nominator parameters. This leads to an interesting and non-trivial relation 
for involved Gauss' functions differing in 0(e) as e — > 0, 



2fl(M; i + £ ; i ) + !fl( £, 1; i ; i). 2 .3-«|ia. 

The second check is given by the relation 



d(d-2)(d-A) 



C(l) 



(30) 



(31) 



where C(l) is the small-g 2 expansion coefficient appearing in eq. (31) of [34]. This coefficient is 
expressible in terms of C(0) by the recurrence relation (37) of [34], 



C(l) 



d- 3 r 



d-2 



d + 4 



C(0) 



(32) 



Now, taking into account the first equality of (29), we reproduce from (31)-(32) the previously 
obtained expression (23) for j(d) in terms of i(d). Thus, the formula (23) is a consequence of 
recurrent relations that hold true for the small-g 2 expansion coefficients of the Feynman inte- 
gral J(m 2 ,q 2 ;d). In turn, these recurrences follow from the differential equation 3 satisfied by the 
function J(m 2 ,q 2 ;d) [34]. 

A full e expansion of the functions i(d) and j(d) can be obtained with the help of eq. (4.16) 
from [37]. Using it we can write 



i(d) 



1 I i £ 27rr(e) e i s e 

I _|_ 3^2^2 ~ ^ 



2-e 



T 2 (e/2) 



l 3 *~ f E$[WHi(y)-I*i + i(*)]} 03) 



where 



Ls 



■(*) = - f 
Jo 



dOlrj- 1 



2 sin ■ 



are the log-sine functions (see [37,56]). Per definition, 

Lsi(0) = -0 and Ls 2 (6») = C\ 2 {6) 



(34) 



(35) 



where Cl 2 (^) is the Clausen's function. This special function appeared readily in Sec. 4, in the 
calculation of the integrals i(d) and j(d) in two dimensions. With the first two terms from the sum, 
we get the e expansions 



lee 2 e 2 



i(d) = o + 7 + V 



Cl 2 (|)+0(e 3 ) 



and 



2 4 8 2V3 V3 

^ = -l-^ + 37! cl2 © + ^ 3 ) 



(36) 
(37) 



3 For a recent review on applications of differential equations in calculations of Feymnan integrals see [41]. 
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In deriving these we took into account that Cl2(27r/3)=2/3-Cl2(7r/3), Cl2(7r)=0, and used the 
relation (23) between i(d) and j{d). The non-trivial 0(e 2 ) term appearing in curly brackets of (33) 
agrees with that of equation (33) in [34]. 

We see that a very special transcendental constant Cl2(7r/3), which is the maximum value of 
the Clausen's integral Cl2(0), appears both in the epsilon expansion of i(d) near d — 4 and in 
the explicit expression of i(d) at d = 2. This is quite natural because the Feynman integrals of 
the type considered here obey certain general relations that connect their values in different space 
dimensions [59,35,39,40]. 

A simple formula relating the values of Feynman integrals associated with the zero-momentum 
"sunrise" diagram (^^) with arbitrary masses on the lines at d — 4 — e and d = 2 — e has been 

found in [59]. In the particular equal-mass case with m 2 = 1, this formula reads, in our notation, 



J(l,0;4-e) 



3tt 2 



(l-e)(2-e) 



J(l,0;2-e)-7r 2 - £ r 2 



(I) 



(38) 



Note that the pole terms of the Laurent expansion of J(4 — e) are contained only in the last, trivial 
term, and the value J(2) directly enters the finite part of J(4 — e) (cf. [39], Sec. 3-4). This feature 
directly maps onto the connection between the functions z(4 — s) and i(2 — e), which can be written 

as 



i(4- 



2-e 



3e 2 



4(1 



■i(2- 



(39) 



It is straightforward to check that the last relation is indeed satisfied by the function i(d) given, 
for instance, by (29). 

Let us return to the epsilon expansions (36)-(37). They can be compared with that quoted in 
(8)-(9). In (36) we see the 0(e 2 ) term of i(d), which is missing in (8). By comparing the 0(e 2 ) 
terms of j(d) we can make the identification 



r 1 dx f 
~ J 1 - x(l - X) \ 



1 + 



ln[x(l - x)} 
[l-x(l-x)} 



8 



3V3 



CI; 



(40) 



An attempt to calculate this integral with the help of Mathematica [60] yields an expression in 
terms of tp'(x) where *p(x) is the psi function, the logarithmic derivative of the Gamma function 
(see e.g. [61]). Thus we get the chain of equalities 



3 r 2 



1113 3 1 



Cl 2 



2V3 V3 



,2' 2' 2' 2' 2' 4J V3 
The rightmost connection can be also read off from (A. 13) of [11]. 



3V3 ' 



(41) 



6. Concluding remarks 



Let us return to the two- loop RG functions defined at the beginning of the Section 2. We see 
that they can be expressed in terms of a single function of d = 4 — e, say 



i(d) 



1 



3 2 Fi( 



£ - l-i-I 

2' '2'4 



-3-i-f 27rr(e) 



P( £ /2) 



(42) 



which contains only one non-trivial Gauss' hypergeometric function. For example, the function 
70 (u) is thus given by 



70 (u) = -4 



4-ri A^ + 2 
3d (iV + 8) 2 



4-(4 + d)i(d) u 2 + 0(u 3 ). 



(43) 



Similar expressions can be constructed also for the remaining RG functions (3(u) and 7^2 (u) . 

Everywhere the dependence on space dimension d is given explicitly in a simple parametric 
form. Obviously, the convenience of such expressions is much better compared to representations 
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in terms of multiple integrals as that in (6)-(7). In principle, this kind of formulas could be used 
in considering some dimensional expansions not tied to d = 4. 

Analytical two-loop results are interesting in view of availability of numerical tables of three- 
loop diagrams of the massive field theory in d dimensions [62]. It could be thought of extending 
explicit calculations to that order. 

Finally, they are also interesting in their own right, as any mathematical results derived in 
closed form. Moreover, in the course of the present work it was possible to find out some interesting 
mathematical relations given in equations (26), (41) and (30). 



A. Explicit results in two dimensions 



In order to derive the values of the functions i(d) and j(d) at d = 2 we have to take e = 2 — a 
and consider the corresponding limits with a — > 0. Thus we write i(d) in (19) as 



i(d) = ii(a) + i 2 (a) 



with 



and 



ii(a) 



2 r l 



(2-a,l; 



3 a 1 

2 ~ 2' 4 



2^1 (2 



a ^ 3 1 

2' '2' 4 



1 „ / a 
12(a) = - 2 F 1 (2--,l 

There exists a simple a — > limit of 12(0:), 



a _ 3 1 

' 2' 4 



,^ 1 2vr 

*2(0) = - + — F 

3 9V3 



(44) 
(45) 
(46) 

(47) 



while each of the hypergeometric functions in i\(a) has to be expanded to 0(a) in order to give the 
finite zi(0). To this end, we found it useful to transform these functions via the linear transformation 
formula ([53], entry 15.3.3) 



2F1 (a, b;c;z) = (l- zf- a - h 2 F l (c-a,c-b; c; z) . 



Thus we get 



*i 



3 X -3/2+a/2 1 

a 



( 1 a 1 a 3 a 1\ / 
2fl V"2 + 2'2" 2 ; 2" 2 ; lJ ~ 2jFi V 



1 a 1.3.1^ 

2 + 2'2' 2' 4/ 



(48) 



• (49) 



Here the both Gauss functions are of the type (11), and hence can be represented in terms of 
simple integrals appearing in (10). We have 



3x-3/2+a/2 1 



2a 



/3x-3/2+a/2 1 

[(l-a)7 1 -/ 2 ] = ( i ) - [{h - h) - ah] 



where 



h = [ dtt~2-f (l _t/4)2-f and 7 2 = / dtt~^(l -t/4)^~^ 
Jo Jo 

At a = 0, the integrals I\ and I 2 both are equal to 

T _ * , V3 
7 °- 3 + ~' 

while the difference between them 1\ — I 2 is of order a: 



h-h 



dt , . , . 1 a . w . a f dt , . . . 1 , „, 9, 

. _(l-t/4)*-T( t -T_i) = __y -^=(l-t/4)Mnf + 0( a 2 ). 



(50) 
(51) 

(52) 
(53) 
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Using Mathematica [60] we get 



7T v/3 , „/l 1 13 3 b 
6 + T 



0(a 2 ). (54) 



f (L L Li i.L\ 

3 2 V2' 2' 2' 2' 2' 4/ 

Combining (50), (54), and (52) we obtain 

2tt1 4 /l 1 1 3 3 1\ 
^--^!-3 + vI 3F2 (2'2'2^'2 ; 4)' < 55 > 

and finally, together with 12(0) from (47), we arrive at i(2) given in equation (24) of the main text. 

We did not find in mathematical tables any summation formula for the generalized hypergeo- 
metric function 3F2 appearing in the last two equations. So we looked for the series representations 
of this function 



V ^ k 1 = IV W [ ( l \ 2k+1 = V ( 2fc - 1 ) !! 1 (561 
^fc!(2fc + l) 2 4 fe (fc!) 2 (2fc + 1) 2 U7 (2fc)!!(2fc + l) 2 4 fc ' v ' 



The middle version of the series appears, indeed, as a partial case of the formula 5.2.13.7 in 
[61] expressed in terms of the function arcsin 2x and its derivative. But, unfortunately, a simple 
numerical check at x = 1/4 revealed the incorrectness of this entry. 
In order to find the right answer, we considered the series expansion 

a(x)-T (V2)fc (57) 
g[X) ~ ^ fc!(2fc + l) 2 4 fc ' (57j 

fe>o v ; 

such that g{\) = 3 F 2 (l/2, 1/2, 1/2; 3/2, 3/2; 1/4) and g(0) = 0. Its derivative with respect to x, 
multiplied by x, is 



^ (l/2) fc x 2 ^ 1 
^W-A.fc! (2 fc + 1) 4* ' 

Differentiating once again we get 



(58) 



, ^ (l/2) fc x 2fc 1 
x (x) = > — 77^ — -r- = = . (59) 

Integrating back leads us to 

X 

xg'(x) — 2 arcsin — (60) 

and 

/2 x xx 

dx— arcsin — = Cl2(2 arcsin — ) + arcsin — lnx (61) 

where Cl2(^) is the Clausen function [56]. In deriving the last equality we used the entry 1.7.4.13 
from [61] and the fact that the integration constants in the both last indefinite integrations vanish. 
Here we observe also that the neighboring entry 1.7.4.12 of [61] involves the series expansion of the 
same type as the rightmost one in (56), giving, along with 1.7.4.11, the corrected version of the 
false formula 5.2.13.7 in the came reference. 

Now, recalling that at x = 1, the function g(x) reproduces the hypergeometric series 3F2 from 
(55), we come to the relation (26) given in the main text. 
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